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It is shown that the previously noted extreme perturbative instability of the longitudinal structure function 
F L (x,Q 2 ) in the small Bjorken-a; region, is a mere artefact of the commonly utilized 'standard' gluon distributions. 
In particular it is demonstrated that using the appropriate - dynamically generated - parton distributions at NLO 
and NNLO, Fl{x, Q 2 ) turns out to be perturbatively rather stable already for Q 2 > O (2 — 3 GeV 2 ). 



A sensitive test of the reliability of perturba- 
tive QCD is provided by studying |1|2|3|4|5] the 
perturbative stability of the longitudinal struc- 
ture function Fl(x, Q 2 ) in the very small Bjorken- 
x region, x < 10~ 3 , at the perturbatively rele- 
vant low values of Q 2 > 0(2 - 3 GeV 2 ). For 
the perturbative-order independent rather flat 
toy model parton distributions in pQ, assumed 
to be relevant at Q 2 ~ 2 GeV 2 , it was shown 
that next-to-next-to-leading order (NNLO) ef- 
fects are quite dramatic at x ;$ 1CV 3 (cf. Fig. 4 
of pQ). To some extent such an enhancement is 
related to the fact, as will be discussed in more 
detail below, that the third-order a 3 contribu- 
tions to the longitudinal coefficient functions be- 
have like xc L ~ — In a; at small x, as compared 
to the small and constant coefficient functions at 
LO and NLO, respectively. It was furthermore 
pointed out, however, that at higher values of Q 2 , 
say Q 2 ~ 30 GeV 2 , where the parton distribu- 
tions are expected to be steeper in the small-a; 
region (cf. eq. (13) of [T]), the NNLO effects are 
reduced considerably. It is well known that dy- 
namically generated parton distributions [5] are 
quite steep in the very small-a; region already at 
rather low Q 2 , and in fact steeper [7] than their 
common 'standard' non-dynamical counterparts. 
Within this latter standard approach, a full NLO 
(2-loop) and NNLO (3-loop) analysis morevover 
confirmed [213] the perturbative fixed-order in- 
stability expectations of [1] in the low Q 2 region. 

It is therefore interesting to study this issue 



concerning the perturbative stability of Fl (x,Q 2 ) 
in the low Q 2 region, Q 2 < 5 GeV 2 , within 
the framework of the dynamical parton model 
|6l7j . For this purpose, following [5], we re- 
peat our previous [8] 'standard' evaluation of the 
NLO and NNLO distributions within the dynam- 
ical approach where the parton distributions at 
Q > 1 GeV are QCD radiatively generated from 
-ya/ence-like (positive) input distributions at an 
optimally determined Qq = \i < 1 GeV (where 
'valence-like' refers to a/ > for all input dis- 
tributions xf(x,fi 2 ) ~ x a f(l — x) bf ). This more 
restrictive ansatz, as compared to the standard 
approach, implies of course less uncertainties [7] 
concerning the behavior of the parton distribu- 
tions in the small— a; region at Q > fi which is 
entirely due to QCD dynamics at x < 10~ 2 . The 
valence-like input distributions at Qq = fl < 1 
are parametrized according to [518) 

xq v (x,Ql)=N qv x a *>{l - x) b '"(l + c g „VS 
+d qv x + e^x 1 ' 5 ), 

xw(x,Ql) = N w x a ™(l-x) b ™(l + c w y/x' + d w x)(l) 

for the valence q v =u v , d v and sea w = q, g den- 
sities, and a vanishing strange sea at Q 2 = Qq, 
s(x, Qq) = s(x, Qq) = 0. All further theoreti- 
cal details relevant for analyzing F^ at NLO and 
NNLO in the MS factorization scheme have been 
presented in [8]. The heavy flavor (dominantly 
charm) contribution to F 2 is taken as given by 
fixed-order NLO perturbation theory 9 10] us- 
ing m c — 1.3 GeV and m& = 4.2 GeV as implied 
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Table 1 

Parameter values of the NNLO and NLO QCD fits with the parameters of the input distributions referring 
to (1) at a common input scale Q\ — fi 2 — 0.5 GeV 2 optimal at both perturbative orders. 
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by optimal fits [7 to recent deep inelastic c- and 
6-production HERA data. Since a NNLO calcu- 
lation of heavy quark production is not yet avail- 
able, we have again used the same NLO 0(a 2 s ) 
result. Finally, we have used for our fit-analyses 
the same deep inelastic HERA-HI, BCDMS and 
NMC data, with the appropriate cuts for Ff as 
in [8] which amounts to a total of 740 data points. 
The required overall normalization factors of the 
data turned out to be 0.98 for HI and BCDMS, 
and 1.0 for NMC. We use here again solely deep 
inelastic scattering data since we are mainly inter- 
ested in the small-x behavior of structure func- 
tions. The resulting parameters of the NLO and 
NNLO fits are summarized in Table 1. The dy- 
namical gluon and sea distributions, evolved to 
some specific values of Q 2 > Q§, are at the NLO 
level very similar to the ones in [7] which were 
obtained from a global analysis including Teva- 
tron Drell-Yan dimuon production and high-Ey 
inclusive jet data as well. Furthermore, the dy- 
namically generated gluon is steeper as x — > 
than the gluon distributions obtained from con- 
ventional 'standard' fits |7l8j (based on some ar- 
bitrarily chosen input scale Qq > 1 GeV 2 , i.e. 
Qq ~ 2 GeV 2 ). On the other hand, the dynam- 
ical sea distribution has a rather similar small-x 
dependence as the 'standard' ones [718] ; this is 
caused by the fact that the valence-like sea input 
in (2) vanishes very slowly as x — > (correspond- 
ing to a small value of dq, a,q ~ 0.05 according 
to Table 1) and thus is similarly increasing with 
decreasing x down to x ~ 0.01 as the sea input 



obtained by a 'standard' fit. Similar remarks hold 
when comparing dynamical and standard distri- 
butions at NNLO. At NNLO the gluon distribu- 
tion xg is flatter as x decreases and, in general, 
falls below the NLO one in the small-cc region, 
typically by 20 - 30% at x ~ 10" 5 and Q 2 < 10 
GeV 2 , whereas the NNLO sea distribution xq is 
about 10 - 20% larger (steeper) than the NLO 
one. 

Now we turn to the perturbative predictions for 
Fl(x, Q 2 ) which can be written as 

2 

x~ l Fi = C L . ns ®q n s+g {Cl, q ®q s + C L . g <g> g) 
+x- 1 F£ (2) 

where £3 in the rif = 3 light quark flavor sec- 
tor denotes the common convolution, q ns stands 
for the usual flavor non-singlet combination and 
9s = J2 q =u,d.s(.q + q) i s tne corresponding flavor- 
singlet quark distribution. Again we use the NLO 
expression [9|10j for F£ also in NNLO due to our 
ignorance of the 0(a 3 s ) NNLO heavy quark cor- 
rections. The perturbative expansion of the coef- 
ficient functions can be written as 

C L Aa s ,x) = Y: (^p-Yc { L %). (3) 

n— 1 ^ ' 

In LO, c« s = fx, eg, s = 0, c« = 241(1 - x) 
and the singlet-quark coefficient function is de- 
composed into the non-singlet and a 'pure sin- 
glet' contribution, c£J = c^ ns + c^. Suf- 
ficiently accurate simplified expressions for the 
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Figure 1. The gluonic contribution F£ to Fl in 
@ with F[ = |a; CL, g ®g in the dynamical (dyn) 
and standard (std) parton approach at NNLO 
and NLO for two representative low values of Q 2 . 
The standard parton distributions utilized in the 
lower panel are taken from [5]. 



NLO and NNLO coefficient functions 



„( 2 ) 



and 



„(3) 



respectively, have been given in [T] 



-L,i 

. It has 



been futhermore noted in [T] that especially for 
Ch.g both the NLO and NNLO contributions are 
rather large over almost the entire x-range. Most 
striking, however, is the behavior of both Cl,q 
and Cl,q at very small values of x: the vanish- 
ingly small LO parts {xc L X \ ~ x 2 ) are negligible 
as compared to the (negative) constant NLO 2- 
loop terms, which in turn are completely over- 
whelmed by the positive NNLO 3-loop singular 

(3) 

corrections xc Li ~ —In a;. This latter singular 
contribution might be indicative for the pertur- 
bative instability at NNLO pQ, as discussed at 
the beginning, but it should be kept in mind that 
a small-x information alone is insufficient for re- 
liable estimates of the convolutions occurring in 
^ when evaluating physical observables. 
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Figure 2. Dynamical parton model NNLO and 
NLO predictions for Fl(x,Q 2 ). 

The gluonic contribution F[ to Fl in ([2]) is 
shown in Fig. [1] at two characteristic low values 
of Q 2 . Although the perturbative instability of 
the subdominant quark contribution F^ as ob- 
tained in a 'standard' fit does not improve for the 
dynamical (sea) quark distributions 5 , it is ev- 
ident from Fig. [1] that the instability disappears 
almost entirely for the dominant dynamical gluon 
contribution already at Q 2 ~ 2 GeV 2 . This im- 
plies that the dynamical predictions for the total 
Fl(x,Q 2 ) become perturbatively stable already 
at the relevant low values of Q 2 > 0(2 - 3 GeV 2 ) 
as shown in Fig. [21 in contrast to the 'standard' 
results in Fig. [3l In the latter case the stability 
has not been fully reached even at Q 2 = 5 GeV 2 
where the NNLO result at x = 10 -5 is more than 
20% larger than the NLO one. A similar dis- 
crepancy prevails for the dynamical predictions 
in Fig. 3 at Q 2 = 2 GeV 2 . This is, however, 
not too surprising since Q 2 = 2 GeV 2 represents 
somehow a borderline value for the leading twist- 
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Figure 3. As in Fig. [2] but for the common stan- 
dard parton distributions as taken from [8]. 

2 contribution to become dominant at small x val- 
ues. This is further corroborated by the observa- 
tion that the dynamical NLO twist-2 fit slightly 
undershoots the HERA data for F 2 at Q 2 ~ 2 
GeV 2 in the small-x region (cf. Fig. 1 of [7J). The 
NLO /NNLO instabilities implied by the standard 
fit results obtained in |2|3j at Q 2 < 5 GeV 2 are 
even more violent than the ones shown in Fig. [3] 
This is mainly due to the negative longitudinal 
cross section (negative Fl(x, Q 2 )) encountered in 
[213] . The perturbative stability in any scenario 
becomes in general better the larger Q 2 , typically 
beyond 5 GeV 2 [IBB] , as shown in Figs. M and 
[3l This is due to the fact that the Q 2 -evolutions 
eventually force any parton distribution to be- 
come sufficiently steep in x. 

To summarize, we have shown that the extreme 
perturbative NNLO/NLO instability of the lon- 
gitudinal structure function Fl at low Q 2 , noted 
in [2-4], is an artefact of the commonly utilized 
'standard' gluon distributions rather than an in- 



dication of a genuine problem of perturbative 
QCD. In fact we have demonstrated that these 
extreme instabilities are reduced considerably al- 
ready at Q 2 — 2 — 3 GeV 2 when utilizing the 
appropriate, dynamically generated, parton dis- 
tributions at NLO and NNLO. These latter par- 
ton distributions have been obtained from a NLO 
and NNLO analysis of F^'" data, employing the 
concepts of the dynamical parton model. It is 
gratifying to notice, once again, the advantage of 
the dynamical parton model approach to pertur- 
bative QCD. 
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